Transient-state control method for a hybrid drive system for vehicles

ABSTRACT

The invention relates to a method of controlling a hybrid drive system having at least one electric machine and at least one diesel engine. The invention corrects the static torque setpoint T eng,s   sp  of the diesel engine and is performed in transient state as follows: constructing a nitrogen oxides (NOx) emission model allowing estimation of the emissions of the diesel engine as a function of torque T eng  of the diesel engine; calculating a dynamic NOx emission setpoint NO x   sp  from the NOx emission model and from the static torque setpoint T eng,s   sp ; determining a dynamic torque setpoint of the diesel engine T eng   sp  from NOx emission setpoint NO x   sp ; determining a dynamic torque setpoint of the electric machine T mot   sp  as a function of the static torque setpoint T pwt   sp  of the drive system and of the dynamic torque setpoint of the diesel engine T eng   sp I and applying the dynamic torque setpoints to the drive system.

CROSS-REFERENCE TO RELATED APPLICATION

Reference is made to French patent application Ser. No. 11/03.497, filed on Nov. 17, 2011, which application is incorporated herein by reference in its entirety.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to an engine control and more particularly to the transient-state control of a hybrid drive system for vehicles, in order to reduce emissions in a hybrid vehicle comprising at least one electric machine and at least one thermal engine for driving the vehicle.

2. Description of the Prior Art

Accounting for emissions in transient state (generated by the change from one steady-state operating point to another steady-state operating point) from an internal-combustion engine is a very important objective. Indeed, emissions from current engines are highly sensitive to operating deviations. Furthermore, future type approval driving cycles will impose more transient phases. In the case of a hybrid vehicle equipped with a diesel engine, the problem is twofold. First, a compromise has to be reached between the NOx emissions (nitrogen oxides, pollutants emitted by the thermal engine) and the fuel consumption for the steady-state operating point, and second, the transient operating phases of the engine, during which NOx emission peaks are observed, have to be managed.

An example illustrating the transient-state emissions problem is given in FIGS. 1 a) to 1 d). This example shows the first acceleration of the extra-urban part of an NEDC type driving cycle (New European Driving Cycle, which is a motor vehicle drive cycle intended to mimic in a reproducible manner the conditions encountered on European roads, mainly used to measure vehicle consumption and emissions) for a hybrid drive system equipped with an exhaust gas recirculation system (EGR). During gear ratio change phases, injection cut-off necessarily leads to a decrease in the richness at the exhaust. Burnt gas availability at the intake is therefore limited when injection is restored upon acceleration. In this case, convergence of the burnt gas fraction at the intake to its setpoint is not instantaneous. Its response time depends on the burnt gas transport time in the EGR circuit. During transient phases, the composition of the gases at the intake (and therefore in the cylinder) does not converge to its setpoint and the NOx emissions exceed the static-state emissions. The amplitude of these emission peaks is high in relation to the expected levels under steady conditions and these peaks represent, according to the type of driving cycle, a not insignificant part of the cumulative emissions.

Over a European driving cycle (NEDC), the transient part does not exceed 15% of the total NOx emissions. On the other hand, for an American driving cycle of FTP type (Federal Test Procedure, a drive cycle used in the USA to measure emissions), the transient part represents more than 40% of the total emissions, as can be seen in Table 1. This table shows a comparison between the measured NOx emissions and the quasi-static NOx emissions for two driving cycles (NEDC and FTP) and the transient part of the total emissions for each cycle.

TABLE 1 DRIVING CYCLE NEDC FTP Exper- Transient Exper- Transient iment Static part iment Static part NOx [mg/km] 109 90 15% 170 100 42% Consumption 4.5 4.3  3% 4.7 4.4  6% [l/100 km]

Accounting for the transient part therefore is an important issue regarding engine control. In this context, thermal engine hybridization represents an option with a high potential. The objective is limiting of thermal engine stress under transient conditions using the electric machine or an alternator starter for compensating for the thermal engine torque setpoint.

Energy supervisors have been developed to control such a hybrid drive system. There are several methods based on optimal control. This type of energy supervisor is presented in the following documents:

-   A. Sciarretta, M. Back & L. Guzzella, “Optimal Control of Parallel     Hybrid Electric Vehicles”, IEEE Transactions on Control Systems     Technology, vol. 12, no. 3, May 2004. -   A. Chasse, A. Sciarretta & J. Chauvin, “Online Optimal Control of a     Parallel Hybrid with Costate Adaptation Rule”, Proceedings of the     IFAC Symposium on Advances in Automotive Control, Munich, Germany,     Jul. 12-14, 2010. -   O. Grondin, L. Thibault, Ph. Moulin, A. Chasse & A. Sciarretta,     “Energy Management Strategy for Diesel Hybrid Electric Vehicle”,     Proceedings of the 7th IEEE Vehicle Power and Propulsion Conference,     Chicago, USA, 6-9 Sep. 2011.

However, although these energy supervisors make it possible to manage the static states of the hybrid drive system, they are not suited for the management of transient states.

The thermal engine is considered to be a purely quasi-static system. This hypothesis is acceptable if the engine torque is considered but it is wrong if the emissions are considered. Thus, the part of the emissions produced during transient operation phases in relation to the total emissions is not insignificant. It is therefore necessary to account for the transition between two steady points selected by the static supervisor.

FIGS. 2 a) and 2 b) show the impact of the torque gradient on the NOx emissions. In order to manage the transient operation phases of the engine, heuristic methods intended to limit the value of the torque gradient required from the thermal engine, while using the electric machine to provide the missing torque and to meet the wheel setpoint, have been developed. This type of heuristic method is illustrated in the following documents:

-   N. Lindenkamp, C.-P. Stöber-Schmidt & P. Eilts, “Strategies for     Reducing NO_(x) and Particulate Matter Emissions in Diesel Hybrid     Electric Vehicles”, SAE Paper n° 2009-01-1305, 2009, -   O. Predelli, F. Bunar, J. Manns, R. Buchwald and A. Sommer, “Laying     Out Diesel-Engine Control Systems in passenger-Car Hybrid Drives”,     Proceedings of the IAV conference on Hybrid Vehicle and Energy     Management, pp. 131-151, Feb. 14-15, 2007.

However, this type of empirical strategy does not allow finding the optimum trajectory and adjustment of the gradient slope depends on the engine, its calibration and the drive mode. In fact, the limits of these gradients have to be calibrated for each cycle.

SUMMARY OF THE INVENTION

The invention allows reduction of the transient emissions of a thermal engine in a hybrid drive system by integrating the physical phenomena involved in their production.

The invention is a method of controlling a hybrid drive system for vehicles in which the torque trajectory that minimizes thermal engine emissions is defined. During transient operation phases, the invention modifies the distribution among the torque setpoint of the electric machine and the torque setpoint of the thermal engine. The principle is to limit the dynamics of the torque setpoint of the thermal engine. The electric machine provides torque compensation to meet the wheel request. Therefore, according to the invention, the NOx emissions of the diesel engine and the hybrid drive system are modelled with a computer and this model is used to define the torque setpoints of the components of the hybrid drive system.

The advantage of this method lies in the fact that the physical phenomena that governs the production of emissions are accounted for in the calculation of the torque trajectory. Using models of the physical phenomena involved ensures implementation of the method and simplifies the adjustment thereof since the method requires calibration of a single adjustment parameter. Furthermore, the control method allows reduction of the transient part of the NOx emissions, while keeping the fuel consumption gains linked with hybridization. The method according to the invention is complementary to conventional (static) strategies and provides uncoupling of the static and transient objectives in the management of hybrid vehicles.

The invention relates to a method of controlling a hybrid drive system having at least one electric machine and at least one diesel engine, wherein a static torque setpoint T_(eng,s) ^(sp) of the diesel engine is acquired, comprising the following:

-   -   constructing a computer-implemented nitrogen oxides (NOx)         emission model allowing estimation of the emissions of the         diesel engine as a function of torque T_(eng) of the diesel         engine;     -   calculating a dynamic NOx emission setpoint NO_(x) ^(sp) from         NOx emission model and from the static torque setpoint T_(eng,s)         ^(sp) of the diesel engine;     -   determining a dynamic torque setpoint of the diesel engine         T_(eng) ^(sp) from NOx emission setpoint NO_(x) ^(sp);     -   determining a dynamic torque setpoint of the electric machine         T_(not) ^(sp) as a function of the dynamic torque setpoint of         the diesel engine T_(eng) ^(sp); and     -   applying the dynamic torque setpoints to the drive system.

In an embodiment the nitrogen oxides emission model depends on the diesel engine speed and on the maximum temperature in the engine cylinders.

Advantageously, the diesel engine is equipped with an exhaust gas recirculation system causing burnt gas to flow into an intake manifold, upstream from a cylinder of the diesel engine. The NOx emission model depends on the burnt gas fraction F₁ ^(est) in the intake manifold of the diesel engine.

Preferably, the NOx emissions are modelled by the formula:

${NO}_{x} = {{\alpha_{1}\left( \frac{N_{e}}{\alpha_{2}} \right)}^{\alpha_{3}}\left( {\alpha_{4}\left( {{\hat{T}}_{cyl} - \alpha_{5}} \right)} \right)^{\alpha_{6}{({1 - {\alpha_{7}F_{1}}})}}}$

with:

-   -   α₁ to α₇ being calibration coefficients of the NOx model,     -   N_(e) being diesel engine speed,     -   {circumflex over (T)}_(cyl) being maximum temperature in the         engine cylinder,     -   F₁ being burnt gas fraction in the intake manifold if the diesel         engine is equipped with an exhaust gas recirculation circuit,         -   this value is zero otherwise, and     -   NO_(x) is the NOx emissions level.

In an embodiment, a NOx emission threshold is defined for the diesel engine. The NOx emission setpoint NO_(x) ^(sp) is calculated as a function of the threshold.

Alternatively, a NOx emission reduction factor ξ can be set with the NOx emission setpoint NO_(x) ^(sp) being calculated as a function of the reduction factor ξ.

Furthermore, determination of the dynamic torque setpoint of the diesel engine can be obtained by inverting the NOx emission model.

Preferably, the static torque setpoint T_(eng,s) ^(sp) of the diesel engine is obtained from a drive system torque setpoint T_(pwt) ^(sp) by an energy supervisor, notably a quasi-static energy supervisor of the hybrid drive system.

The invention also relates to a computer program product downloadable from a communication network and/or recorded on a computer readable medium and/or controller executable. This computer program product comprises program code instructions for implementing the method according to the invention when the program is executed on a controller.

The invention further relates to a hybrid drive system for a vehicle having at least one electric machine and at least one diesel engine. The drive system is controlled by the control method as defined above. Control of the drive system in transient state can be in cascade with respect to a quasi-static energy supervisor.

The invention also relates to a vehicle, notably a motor vehicle, comprising a hybrid drive system according to the invention.

BRIEF DESCRIPTION OF THE DRAWINGS

Other features and advantages of the method according to the invention will be clear from reading the description hereafter of embodiments given by way of non limitative example, with reference to the accompanying figures wherein:

FIGS. 1 a) to 1 d) illustrate a comparison between the measured NOx emissions and the quasi-static NOx emissions for the first acceleration of the extra-urban part of the NEDC driving cycle, in terms of vehicle speed (1 a)), burnt gas fraction in the intake manifold (1 b)), engine speed (1 c)) and NOx emissions (1 d));

FIGS. 2 a) and 2 b) illustrate the influence of the engine torque gradient (dT) on the NOx emissions in transient state;

FIG. 3 illustrates the dynamic setpoint of the thermal engine torque and of the nitrogen oxides emissions;

FIG. 4 is a flowsheet of the dynamic strategy;

FIGS. 5 a) and 5 b) are flowsheets for the calculation of a torque setpoint trajectory for the thermal engine. FIG. 5 a) shows the NOx emission model and FIG. 5 b) illustrates a specific embodiment;

FIG. 6 illustrates the principle of the NOx setpoint calculation;

FIG. 7 illustrates the maximum temperature map of the cylinder;

FIGS. 8 a) to 8 f) show simulation results for the method of controlling a hybrid drive system according to the invention, in terms of thermal engine torque (8 a)), electric motor torque (8 b)), burnt gas fraction (8 c)), cylinder temperature (8 d)), and NOx emissions (8 e)) and (80);

FIG. 9 shows the evolution of the fuel consumption as a function of the NOx emissions for various reduction factor values; and

FIG. 10 illustrates the part of the transient NOx emissions in an FTP cycle as a function of the reduction factor.

DETAILED DESCRIPTION OF THE INVENTION

The method according to the invention allows reduction of NOx emissions for a hybrid drive system. It allows management of transient phases between two static operating points using a NOx emission model that can predict NOx emissions. This model allows determination of the NOx emissions of the hybrid drive system as a function of the torque of the internal-combustion engine.

According to the invention, the method allows controlling a hybrid drive system of a vehicle, notably a motor vehicle, having at least one electric machine and at least one diesel engine. For the transient states, the following stages are carried out:

-   -   acquiring a static torque setpoint T_(eng,s) ^(sp) of the         thermal engine;     -   constructing a NOx emission model allowing estimation of the NOx         emissions of the diesel engine as a function of the torque         T_(eng) of the diesel engine;     -   calculating a dynamic NOx emission setpoint NO_(x) ^(sp) from         the NOx emission model determining the NOx emissions of the         diesel engine as a function of the torque setpoint T_(eng,s)         ^(sp);     -   deducing a dynamic torque setpoint T_(eng,s) ^(sp) of the diesel         engine     -   determining a dynamic torque setpoint T_(mot) ^(sp) of the         electric machine; and     -   applying the dynamic torque setpoints to the drive system.

Nomenclature

The following notations are used in the description hereafter:

-   -   ξ is a NOx peak reduction factor [in %]′     -   β is a NOx setpoint initialization threshold [in ppm]     -   α_(i) are calibration parameters of the NOx emissions model         [unitless]′     -   F₁ ^(sp) is the burnt gas mass fraction setpoint in the intake         manifold [in %];     -   F₁ ^(est) is the estimation of the burnt gas mass fraction in         the intake manifold [in %];     -   N_(e) is the thermal engine speed [in rpm]′     -   T_(pwt) ^(sp) is the static torque setpoint of the drive system         [in Nm];     -   T_(eng,s) ^(sp) is the thermal engine static torque setpoint of         the thermal engine to be reached as a function of the torque         setpoint of the drive system [in Nm];     -   T_(mot,s) ^(sp) is the electric machine static torque setpoint         to be reached as a function of the torque setpoint of the drive         system [in Nm];     -   T_(eng) ^(sp) is the dynamic thermal engine torque setpoint [in         Nm];     -   T_(mot) ^(sp) is the dynamic electric machine torque setpoint         [in Nm];     -   {circumflex over (T)}_(cyl) ^(sp) is the maximum cylinder         temperature for the static torque setpoint [in K];     -   {circumflex over (T)}_(cyl) ^(dyn) is the maximum cylinder         temperature minimizing NOx emissions [in K];     -   ΔNO_(x) is the NOx peak amplitude with respect to the target         steady point [in ppm]     -   εNO_(x) is the NOx emissions error [in ppm];     -   NO_(x) is the maximum value of the NOx emission peak [in ppm];     -   NO_(x) ^(s) is the value of the NOx emissions on a steady point         [in ppm];     -   NO_(x)(t) is the NOx emissions generated by the nominal torque         transient [in ppm];     -   NO_(x) ^(sp) is the NOx emissions setpoint [in ppm].

The principle of the method calculates, during transient operation phases, a thermal engine torque trajectory that minimizes the NOx emissions. The principle of such a method is shown in FIG. 3. It determines, when the torque setpoint ranges between a point A and a point B, if it leads to transient surplus NOx emissions. The curves in thick line show the variation of the thermal engine torque and the variation of the NOx emissions obtained with the method according to the invention and the dotted lines show an instant thermal engine torque variation.

In a preferred embodiment, the drive system is equipped with an exhaust gas recirculation circuit (EGR). In this case, a high engine torque variation generally causes great deviations in the air loop (air intake and exhaust gas recirculation). The thermodynamic conditions in the combustion chamber then deviate from the steady-state conditions and the NOx emissions increase. It is therefore desirable to calculate a new trajectory for the thermal engine torque. The method of determining this new torque setpoint uses a NOx emission model allowing determination of the trajectory that provides minimum emissions between two static torque points (A and B). The transition from point A to point B is no longer instantaneous, it follows a trajectory leading to point B′. During this phase, the thermal engine torque is compensated by the electric machine so that the wheel torque request is met. The method thus compensates the thermal engine torque so as to limit demands leading to NOx emission deviations.

PRESENTATION OF AN EMBODIMENT

The embodiment described below relates to a diesel engine equipped with a burnt gas recirculation circuit. Furthermore, in this embodiment, an energy supervisor, of static or quasi-static type, allows determination of the torque setpoints of the hybrid drive system in steady state.

Stage 1)—Acquisition of a Torque Setpoint

The goal of the invention is to reduce the transient part of the NOx emissions. What is referred to as transient state is the shift from a static torque setpoint to a second static torque setpoint.

As illustrated in FIG. 4, the static torque setpoints to be applied to the thermal engine T_(eng,s) ^(sp) and to the electric machine T_(mot,s) ^(sp) are determined from the torque setpoint of the drive system T_(pwt) ^(sp). These static torque setpoints are determined by an energy supervisor, in combination to which the control method according to the invention is implemented. The purpose of this energy supervisor, notably of quasi-static type, is to select the operating points and therefore the torque setpoints of the thermal engine and of the electric machine. The invention is complementary to an ECMS type quasi-static energy supervisor. However, it is totally independent of it and could work downstream from any energy supervisor. The torque setpoint of the drive system T_(pwt) ^(sp) corresponds to the wheel torque request notably required by the driver of the vehicle.

Stage 2)—Construction of the NOx Emission Model

The dynamic NOx emission setpoint is calculated from a NOx emission model determining the NOx emissions of the diesel engine. What is referred to as a NOx emission model is a model that can predict the NOx emissions of the diesel engine as a function of the torque of the diesel engine of the hybrid drive system. In order to be precise and representative of the engine, this model has to be calibrated to the diesel engine of the drive system used, and it has to depend on the operating conditions thereof, for example the thermal engine speed N_(e). Besides, the adjective “dynamic” makes it clear that it is a setpoint that applies to a transient state, that is between two steady-state operating points, evolving as a function of time.

This model has to be invertible so as to allow to determine thereafter the dynamic torque setpoint from the dynamic NOx emission setpoint. A semi-empirical NOx emission model can be used.

In this embodiment, the torque trajectory is calculated in cascade with respect to the static calculation as shown in FIG. 4. The static torque setpoint determined by means of the static optimization strategy, T_(eng,s) ^(sp), is corrected and the dynamic torque setpoint is denoted by T_(eng) ^(sp). The relations connecting the torque of the thermal engine, of the electric machine and the torque setpoint of the hybrid drive system remain unchanged:

T _(pwt) ^(sp) =T _(eng) ^(sp) +T _(mot) ^(sp)  (1)

T _(pwt,s) ^(sp) =T _(eng,s) ^(sp) +T _(mot,s) ^(sp)  (2)

The NOx emissions are calculated by the semi-empirical model as follows:

NO_(x)=φ(N _(e) ,{circumflex over (T)} _(cy) ,F ₁  (3)

It is a correlation that connects the emissions to the engine speed N_(e), the burnt gas fraction at the intake F₁ and the maximum temperature in the cylinder {circumflex over (T)}_(cyl). These quantities have been selected because they are involved in the NOx production in the combustion chamber of a diesel engine. This model is described in detail in the following document:

-   J.-C. Schmitt, M. Fremovici, O. Grondin & F. Le Berr, “Compression     Ignition Engine Model Supporting Powertrain Development, Proceedings     of the IFAC workshop on engine and powertrain control, simulation     and modeling”, Rueil-Malmaison, France, 30 Nov.-2 Dec. 2009.

Function φ is given by:

$\begin{matrix} {{NO}_{x} = {{\alpha_{1}\left( \frac{Ne}{\alpha_{2}} \right)}^{\alpha_{3}}\left( {\alpha_{4}\left( {\hat{T} - \alpha_{5}} \right)} \right)^{\alpha_{6}{({1 - {\alpha_{7}F}})}}}} & (4) \end{matrix}$

Coefficients α_(i) are calibration parameters of the semi-physical NOx model. The values of these coefficients are calculated from experimental data, for example from tests on the engine test bench, and they are valid for an engine. Generally, the data obtained from the calibration of the thermal engine are sufficient for calibration of this type of model.

Stage 3)—Calculation of the NOx Emission Setpoint

FIG. 5 a) shows the principle for calculating the dynamic NOx emission setpoint of the thermal engine, NO_(x) ^(sp), in transient state in a general context. The method uses a NOx emission model determining the NOx emissions of the diesel engine.

FIG. 5 b) shows the principle of this calculation in the particular case of the embodiment described. The model of the drive system is based on a model representative of the NOx emissions and on a model representative of the engine air loop dynamics. Air loop (10) is understood to be the entire intake circuit and exhaust gas recirculation circuit. The complexity level of these two models is in keeping with an integration into an energy optimization strategy that can be implemented in real time in an electronic calculator. It is a semi-empirical model (or “grey-box” model) for NOx emissions and a first-order filter including a pure delay representative of the dynamics of the burnt gas fraction at the engine intake.

These NOx emission and burnt gas dynamics models depend on the engine speed N_(e) and on the burnt gas fraction setpoint in the intake manifold F₁ ^(sp), upstream from a cylinder of the diesel engine, which can be given by a map depending on engine speed N_(e) and on torque setpoint T_(eng,s) ^(sp).

FIG. 5 b) shows the two main stages in the calculation of the torque trajectory: determination of the NOx emissions discharged under transient operating conditions and calculation of a new trajectory of the engine torque setpoint from a reduced emission setpoint in relation to the nominal case. The NOx model is first supplied with the torque setpoint T_(eng,s) ^(sp) resulting from the static optimization. The model, depending on the current setpoint and not on that of the previous calculation step, determines the NOx peak amplitude that would be generated if it was not corrected. This value serves as the basis for calculating the NOx setpoint in the calculation of the new torque trajectory. To sum up, the NOx model and the air loop dynamics model ensure filtering of the torque setpoint that pertains to the dynamics specific to the system so as to limit the emissions peak.

For this model, a NOx emission attenuation rate is defined by an adjustment parameter ξ corresponding to the amplitude attenuation percentage of the initial NOx peak. It is thus possible to adjust the NOx reduction to the value of parameter ξ. For example, in FIG. 6, the transient between two operating points can lead to exceeding of the steady-state emission level denoted by NO_(x) ^(s). The amplitude of the exceeding the steady-state emission level is denoted by NOx. It is thus possible to adapt the control according to outside conditions by adjusting parameter ξ. For example, in the case of an urban zone where traffic is reserved for vehicles with a low emissions level, the parameter can be adapted in order to meet the objective.

In case of significantly exceeding the steady-state emission level, an adjustment parameter allows defining an acceptable threshold crossing that will serve as the setpoint:

$\begin{matrix} {{NO}_{x}^{sp} = {{NO}_{x}^{s} + {\Delta \; {{NO}_{x}\left( {1 - \frac{\xi}{100}} \right)}}}} & (5) \end{matrix}$

where ΔNO_(x) is the difference between the steady-state emissions level and the maximum transient-state peak value:

ΔNO_(x)= NO _(x)−NO_(x) ^(s)  (6)

where NO_(x) ^(s) is the NOx emission level under steady state conditions. Parameter ξ allows to adjust the intensity of the transient-state correction. If parameter ξ is zero, then the NOx limitation strategy is inactive. Otherwise, an acceptable value of exceeding the steady-stated emission level is defined.

The NOx emission model, for example the model of Equation (4), is used for calculation of:

-   -   i. The static value of the steady operating point NO_(x) ^(s),     -   ii. The instant value of the emissions generated by the nominal         torque transient NO_(x)(t), which allows determination of the         dynamic NOx emission setpoint NO_(x) ^(sp), and     -   iii. The instant value of the emissions generated by the torque         transient obtained from the limitation strategy NO_(x)         ^(dyn)(t). This value is not useful for the practice of the         method, but it allows having knowledge of the value of the NOx         emissions produced during the transient state.

i. Calculation of the Static Value of the Steady Operating Point NO_(x) ^(s)

The static value of the NOx emissions on the steady-state operating point NO_(x) ^(s) is given by:

NO_(x) ^(s)=φ(N _(e) ,{circumflex over (T)} _(cyl) ^(sp) ,F ₁ ^(sp))  (7)

where {circumflex over (T)}_(cyl) ^(sp) is the maximum temperature in the cylinder calculated from a map (FIG. 7) depending on engine speed N_(e) and on the torque value on the steady-state operating point:

{circumflex over (T)} _(cyl) ^(sp)=ψ(N _(e) ,T _(eng,s) ^(sp))  (8)

and F₁ ^(sp) is the burnt gas fraction setpoint in the intake manifold on the steady-state operating point. It is given by a map depending on the engine speed and on the torque setpoint:

F ₁ ^(sp) =f(N _(e) ,T _(eng,s) ^(sp))  (9)

By combining Equations (7) to (9), the value NO_(x) ^(s) of the NOx emissions on the steady-state operating point as a function of the static torque setpoint of the diesel engine T_(eng,s) ^(sp) and of the diesel engine speed N_(e) is obtained.

ii. Calculation of the Instant Value of the Emissions Generated by the Nominal Torque Transient NO_(x)(t)

The instant value of the emissions generated by the nominal torque transient NO_(x)(t) is given by:

NO_(x)(t)=φ(N _(e) ,{circumflex over (T)} _(cyl) ^(sp) ,F ₁ ^(est)(t))  (10)

where {circumflex over (T)}_(eyl) ^(sp) is the maximum temperature in the cylinder (same value as for the calculation of NO_(x) ^(s), Equation (8)) and F₁ ^(est)(t) is the estimated burnt gas fraction in the intake manifold. It involves the burnt gas dynamics at the intake. The gas transport dynamics model used is a first-order filter associated with a pure delay such as:

$\begin{matrix} {{{{\tau (t)}\frac{F_{1}^{{est}{(t)}}}{t}} + {F_{1}^{est}(t)}} = {F_{1}^{sp}\left( {t - {t_{r}(t)}} \right)}} & (11) \end{matrix}$

where the filter time constant r and the pure delay t_(r) are parametrized as a function of the engine speed:

$\begin{matrix} {{t_{r}(t)} = \left\{ {\begin{matrix} \frac{k_{r}}{N_{e}(t)} & {{{if}\mspace{14mu} {T_{eng}(t)}} \succ 0} \\ 0 & {otherwise} \end{matrix}{and}} \right.} & (12) \\ {{\tau (t)} = \frac{k_{\tau}}{N_{e}(t)}} & (13) \end{matrix}$

Coefficients k_(r) and k_(τ) allow the burnt gas dynamics to be adjusted and they are two parameters for adjusting the transient emissions limitation strategy. They can be determined experimentally, for example during tests on the engine test bench. The instant value of the emissions generated by the nominal torque transient NO_(x)(t) is involved in the calculation of the amplitude of the NOx peak in relation to the steady-state point ΔNO_(x) described by Equation (4). The value of NO _(x) is calculated from the maximum amplitude of the transient emission peak and it is fixed as long as error ENO_(x) is above a set threshold β. The quantities εNO_(x) represents the error between the emissions generated by the nominal torque transient NO_(x)(t) and the value of the steady point NO_(x) ^(s). Once the system is in the steady state, it is necessary to reinitialize NO _(x) to the static value NO_(x) ^(s). This condition guarantees the independence of two consecutive transients. Finally, the maximum value of the NOx peak is written as follows:

$\begin{matrix} {{\overset{\_}{NO}}_{x} = \left\{ {\begin{matrix} {\max \left\{ {{NO}_{x}(t)} \right\}} & {{{if}\mspace{14mu} ɛ\; {NO}_{x}} \succ \beta} \\ {NO}_{x}^{s} & {otherwise} \end{matrix}{with}} \right.} & (14) \\ {{ɛ\; {NO}_{x}} = {{{NO}_{x}(t)} - {NO}_{x}^{s}}} & (15) \end{matrix}$

By integrating into Equation (6) this value of NO _(x) and value NO_(x) ^(s) of the NOx emissions on the steady operating point determined in the previous stage, ΔNO_(x) is calculated, then the dynamic NOx emission setpoint NO_(x) ^(sp) is determined by applying Equation (5).

iii. Calculation of the Instant Value of the Emissions Generated by the Torque Transient Obtained from the Limitation Strategy NO_(x) ^(dyn)(t)

The instant value of the emissions generated by the torque transient obtained from the limitation strategy NO_(x) ^(dyn)(t) is given by:

NO_(x) ^(dyn)(t)=φ(N _(e) ,{circumflex over (T)} _(cyl) ^(dyn) ,F ₁ ^(est)(t))  (16)

This quantity is not used in the method. However, it provides information on the value of the emissions produced under transient conditions. It is in principle always below the estimated value of the emissions generated by the nominal torque transient. In Equation (16), F₁ ^(est) is the estimated intake composition (same value as for the calculation of NO_(x)(t), Equation (11)) and {circumflex over (T)}_(cyl) ^(dyn) is the maximum temperature in the cylinder calculated from the corrected torque setpoint T_(eng) ^(sp):

{circumflex over (T)} _(cyl) ^(dyn)=ψ(N _(e) ,T _(eng) ^(sp))  (17)

Stage 4)—Determination of the Diesel Engine Dynamic Torque Setpoint

Once the setpoint NO_(x) ^(sp) is calculated, the NOx emission model is inverted to allow calculation of a new dynamic torque setpoint T_(eng) ^(sp). The NOx has two purposes in the strategy. First of all, it allows calculation of the response of the system to a torque setpoint resulting from the static optimization strategy and determining therefrom a NOx setpoint. Then, it is inverted and allows calculation of the torque trajectory meeting this setpoint.

In the embodiment described, the previous equations, notably (6) and (8), can be inverted. The corrected torque trajectory is obtained by inverting the maximum cylinder temperature model described by Equation (8):

T _(eng) ^(sp)=ψ⁻¹(N _(e) ,{circumflex over (T)} _(cyl) ^(dyn))  (18)

Temperature {circumflex over (T)}_(cyl) ^(dyn) corresponds to the maximum temperature trajectory in the combustion chamber such that the emissions meet setpoint NO_(x) ^(sp). Temperature {circumflex over (T)}_(cyl) ^(dyn) is calculated by inverting the emission model given by Equation (5) and from setpoint NO_(x) ^(sp) and the estimation of the gas composition at the intake:

{circumflex over (T)} _(cyl) ^(dyn)=φ⁻¹(N _(e) ,NO _(x) ^(sp) ,F ₁ ^(est))  (19)

Function φ⁻¹ is written as follows:

$\begin{matrix} {{\hat{T}}_{cyl}^{dyn} = {\alpha_{5} + {\frac{1}{\alpha_{4}}\left( {{NO}_{x}^{sp}\frac{1}{\alpha_{1}}\left( \frac{\alpha_{2}}{N_{e}} \right)^{\alpha_{3}}} \right)^{\frac{1}{\alpha_{6}{({1 - {\alpha_{7}F_{1}^{est}}})}}}}}} & (20) \end{matrix}$

By combining Equations (18) and (20), the dynamic torque setpoint of the diesel engine T_(eng) ^(sp) is obtained with this setpoint depending on engine speed N_(e), the estimated burnt gas fraction in the intake manifold F₁ ^(est) and the dynamic NOx emission setpoint NO_(x) ^(sp).

Stage 5)—Determination of the Dynamic Torque Setpoint of the Electric Machine

The method according to the invention provides thermal engine torque compensation under transient conditions by at least one electric machine so that the wheel torque request is met. It is therefore necessary to determine the dynamic torque setpoint of the electric machine. This torque setpoint can be obtained by inverting Equation (1):

T _(mot) ^(sp) =T _(pwt) ^(sp) −T _(eng) ^(sp)  (21)

Stage 6)—Application of the Dynamic Torque Setpoints

The invention allows determination of the torque setpoints of the drive of the hybrid system. By applying these setpoints to the thermal engine and to the electric machine, an emissions decrease is obtained and the fuel consumption can also be limited.

Segmentation of the static and transient optimization layers allows uncoupling into two sampling periods which are high frequency for transient optimization and low frequency for the static part. The dynamic optimization period is suited to the physical phenomena involved in the engine, in the present case the production of emissions.

The parameters of the method and their orders of magnitude are summed up in Table 2.

TABLE 2 List of the adjustment parameters of the strategy Symbol Description Order of magnitude ζ NOx peak reduction factor [0; 100] α_(i) 7 calibration coefficients of Depends on the the NOx emissions model engine considered k_(r) Pure delay adjustment coefficient [200; 1000] k_(τ) Adjustment coefficient of the [500; 2000] first-order filter time constant β NOx setpoint initialization [1; 10] threshold F₁ ^(sp) = f(N_(e), T) Map of the burnt gas fraction Constructed from setpoint at the intake engine calibration data {circumflex over (T)}_(cyl) = ψ(N_(e), T) Map of the maximum cylinder Constructed from temperature calibration data

Variants

The invention has been described above within the context of a preferred embodiment. However, the invention can also relate to different variants that can be combined.

Stage 1)—Acquisition of a Torque Setpoint

The static torque setpoint can directly result from a request from the car driver and not from an energy supervisor of static or quasi-static type, notably in an embodiment where the hybrid drive system is provided with no energy supervisor.

Stage 2)—Construction of the NOx Emission Model

If the thermal engine is not equipped with an exhaust gas recirculation circuit, an invertible NOx emission model depending only on the engine speed and the cylinder temperature is defined as NO_(x)(t)=φ(N_(e),{circumflex over (T)}_(cyl),0)=φ′(N_(e),{circumflex over (T)}_(cyl)). The cylinder temperature being a function of the thermal engine torque setpoint ({circumflex over (T)}_(cyl) ^(dyn)=ψ(N_(e),T_(eng) ^(sp)) remains valid. As in the preferred embodiment, this NOx emission model is a semi-empirical model obtained by correlation with experimental tests on engine test benches. For example, Equation (4) can be rewritten by assuming that the burnt gas fraction in the intake manifold is zero. Thus:

${NO}_{x} = {{\alpha_{1}\left( \frac{Ne}{\alpha_{2}} \right)}^{\alpha_{3}}\left( {\alpha_{4}\left( {{\hat{T}}_{cyl} - \alpha_{5}} \right)} \right)^{\alpha_{6}}}$

Stage 3)—Calculation of the NOx Emission Setpoint

Instead of using adjustment parameter ξ to determine a NOx emission attenuation rate, it is possible to define a NOx emission threshold S that should not be exceeded. The threshold S can be fixed to whatever the static state is to be reached, or it can depend on the target steady operating point. The following equation therefore has to be satisfied: NO_(x) ^(sp)≦S. This equation is introduced in stage ii) in place of Equation (3) of the embodiment described. This threshold can be set according to the legislation in force concerning emissions.

Stage 4)—Determination of the diesel engine dynamic torque setpoint When the diesel engine is not equipped with an exhaust gas recirculation circuit, the NOx emission model is inverted: NO_(x)(t)=φ′(N_(e),{circumflex over (T)}_(eyl)) so as to deduce the diesel engine dynamic torque setpoint: {circumflex over (T)}_(cyl) ^(dyn)=φ′⁻¹(N_(e),NO_(x) ^(sp)) then T_(eng) ^(sp)=ψ⁻¹(N_(e),{circumflex over (T)}_(cyl) ^(dyn)).

Stage 5)—Determination of the Dynamic Torque Setpoint of the Electric Machine

Other types of compensation taking account of the capacities (operating range) of the electric machine can be used in the control method according to the invention. For example, if the dynamic torque setpoint of the electric machine T_(mot) ^(sp) exceeds the maximum torque allowable by the electric machine, the torque setpoint of the electric machine is set at the maximum torque, then the dynamic torque setpoint of the thermal engine T_(eng) ^(sp) is recalculated by adapting it to the dynamic torque setpoint change of the electric machine. The best use of the components of the hybrid drive system is thus made with a non-optimal emissions reduction.

Application Example

The results presented in this part are obtained from simulations of a diesel hybrid vehicle with a parallel architecture. The diesel engine is equipped with a burnt gas recirculation system. The thermal engine model is a physical model that can predict emissions and, in particular, NOx emissions. Such a model is for example described in the following document:

-   R. Lebas, G. Mauvio, F. Le Berr & A. Albrecht, “A Phenomenological     Approach to Model Diesel Engine Combustion and In-Cylinder Pollutant     Emissions Adapted to Control Strategy, Proceedings of the IFAC     workshop on engine and powertrain control, simulation and modeling”,     Rueil-Malmaison, France, 30 Nov.-2 Dec. 2009.

This model has been calibrated from tests carried out on a Euro 6 diesel engine. The strategy of emissions limitation under transient conditions has been coupled with a static energy supervisor whose principle is based on optimal control. It is presented in the aforementioned document:

-   O. Grondin, L. Thibault, Ph. Moulin, A. Chasse & A. Sciarretta,     “Energy Management Strategy for Diesel Hybrid Electric Vehicle”,     Proceedings of the 7^(th) IEEE Vehicle Power and Propulsion     Conference, Chicago, USA, 6-9 September 2011.

FIGS. 8 a) to 8 f) show a portion of an FTP driving cycle where the control method is active. It is first noted that the steady operating points are not modified. On the other hand, in transient state, the torque of the thermal engine is limited (FIG. 8 a)). The limitation operates in phases where there are great differences between the setpoint and the estimation of the burnt gas fraction at the intake (FIG. 8 c)). In this case, the internal NOx model anticipates a significant exceeding of the expected level in steady state (FIG. 8 e)). The torque distribution changes and the electric machine allowed to limit the thermal engine dynamics (FIG. 8 b)). The emission peaks that appear during these phases (FIG. 8 f)) are greatly reduced. It is observed that the emissions peak reduction results from the temperature limitation in the cylinder (FIG. 8 d)). Indeed, this quantity is extremely important in the kinetic NOx production mechanisms. These results show that taking account, in an energy supervisor, of physical quantities that govern the production of NOx emissions leads to good results.

FIG. 9 shows the hybrid vehicle consumption and the NOx emissions obtained in an FTP driving cycle and for an adjustment parameter value scan for these simulations. The results show that it is possible to find values of this parameter allowing decreasing the value of the NOx emissions without modifying the consumption of the vehicle. This is the case when ξ ranges between 10 and 40%. For higher values, the emissions continue to decrease but the consumption can increase substantially. In this case, parameter ξ allows adjustment of the compromise between the NOx emissions and the fuel consumption. In any case, the part of the transient NOx emissions in an FTP cycle decreases as parameter ξ increases as shown in FIG. 10.

Furthermore, the method according to the invention affords the following advantages:

-   -   it makes it possible to modify online the emissions level during         transient phases, use of the physical phenomena models involved         which are a semi-physical model (of correlation type) for NOx         emissions prediction and an air loop dynamics model limiting the         number of adjustment parameters, and     -   this method is suited for a hybrid vehicle equipped with a         compression ignition engine but it could however be suited for a         spark-ignition engine. 

1-12. (canceled)
 13. A computer-implemented method of controlling a hybrid drive system including at least one electric machine and at least one diesel engine, wherein a static torque setpoint T_(eng,s) ^(sp) of the diesel engine is acquired, comprising: constructing a nitrogen oxide (NOx) emission model allowing estimation of the emissions of the diesel engine as a function of a torque T_(eng) of the diesel engine; calculating a dynamic NOx emission setpoint NO_(x) ^(sp) from the NOx emission model and from the static torque setpoint T_(eng,s) ^(sp) of the diesel engine; determining a dynamic torque setpoint of the diesel engine T_(eng) ^(sp) from NOx emission setpoint NO_(x) ^(sp); determining a dynamic torque setpoint of the electric machine T_(mot) ^(sp) as a function of the dynamic torque setpoint of the diesel engine T_(eng) ^(sp); and applying the dynamic torque setpoints to the drive system.
 14. A method as claimed in claim 13, wherein: the nitrogen oxide emission model depends on speed of the diesel engine and on a maximum temperature in the engine cylinder.
 15. A method as claimed in claim 13, wherein: the diesel engine is equipped with an exhaust gas recirculation system providing burnt gas to an intake manifold, upstream from a cylinder of the diesel engine and NOx emission model depends on a burnt gas fraction F₁ ^(est) in an intake manifold of the diesel engine.
 16. A method as claimed in claim 14, wherein: the diesel engine is equipped with an exhaust gas recirculation system providing burnt gas to an intake manifold, upstream from a cylinder of the diesel engine and NOx emission model depends on a burnt gas fraction F₁ ^(est) in an intake manifold of the diesel engine.
 17. A method as claimed in claim 13, wherein: the NOx emissions are modelled by the formula: ${NO}_{x} = {{\alpha_{1}\left( \frac{N_{e}}{\alpha_{2}} \right)}^{\alpha_{3}}\left( {\alpha_{4}\left( {{\hat{T}}_{cyl} - \alpha_{5}} \right)} \right)^{\alpha_{6}{({1 - {\alpha_{7}F_{1}}})}}}$ with: α₁ to α₇ being calibration coefficients of the NOx model; N_(e) being diesel engine speed; {circumflex over (T)}_(eyl) being a maximum temperature in an engine cylinder; F₁ being a burnt gas fraction in the intake manifold of a diesel engine equipped with an exhaust gas recirculation circuit and with F₁ being zero otherwise; and NO_(x) is a NOx emissions level.
 18. A method as claimed in claim 14, wherein: the NOx emissions are modelled by the formula: ${NO}_{x} = {{\alpha_{1}\left( \frac{N_{e}}{\alpha_{2}} \right)}^{\alpha_{3}}\left( {\alpha_{4}\left( {{\hat{T}}_{cyl} - \alpha_{5}} \right)} \right)^{\alpha_{6}{({1 - {\alpha_{7}F_{1}}})}}}$ with: α₁ to α₇ being calibration coefficients of the NOx model; N_(e) being diesel engine speed; {circumflex over (T)}_(cyl) being a maximum temperature in an engine cylinder; F₁ being a burnt gas fraction in the intake manifold of a diesel engine equipped with an exhaust gas recirculation circuit and with F₁ being zero otherwise; and NO_(x) is a NOx emissions level.
 19. A method as claimed in claim 15, wherein: the NOx emissions are modelled by the formula: ${NO}_{x} = {{\alpha_{1}\left( \frac{N_{e}}{\alpha_{2}} \right)}^{\alpha_{3}}\left( {\alpha_{4}\left( {{\hat{T}}_{cyl} - \alpha_{5}} \right)} \right)^{\alpha_{6}{({1 - {\alpha_{7}F_{1}}})}}}$ with: α₁ to α₇ being calibration coefficients of the NOx model; N_(e) being diesel engine speed; {circumflex over (T)}_(cyl) being a maximum temperature in an engine cylinder; F₁ being a burnt gas fraction in the intake manifold of a diesel engine equipped with an exhaust gas recirculation circuit and with F₁ being zero otherwise; and NO_(x) is a NOx emissions level.
 20. A method as claimed in claim 13, wherein: a NOx emission threshold is defined for the diesel engine and the NOx emission setpoint NO_(x) ^(sp) is calculated as a function of the threshold.
 21. A method as claimed in claim 14, wherein: a NOx emission threshold is defined for the diesel engine and the NOx emission setpoint NO_(x) ^(sp) is calculated as a function of the threshold.
 22. A method as claimed in claim 15, wherein: a NOx emission threshold is defined for the diesel engine and the NOx emission setpoint NO_(x) ^(sp) is calculated as a function of the threshold.
 23. A method as claimed in claim 17, wherein: a NOx emission threshold is defined for the diesel engine and the NOx emission setpoint NO_(x) ^(sp) is calculated as a function of the threshold.
 24. A method as claimed in claim 13, wherein: a NOx emission reduction factor ξ is set and the NOx emission setpoint NO_(x) ^(sp) is calculated as a function of the reduction factor ξ.
 25. A method as claimed in claim 14, wherein: a NOx emission reduction factor ξ is set and the NOx emission setpoint NO_(x) ^(sp) is calculated as a function of the reduction factor ξ.
 26. A method as claimed in claim 15, wherein: a NOx emission reduction factor ξ is set and the NOx emission setpoint NO_(x) ^(sp) is calculated as a function of the reduction factor ξ.
 27. A method as claimed in claim 17, wherein: a NOx emission reduction factor ξ is set and the NOx emission setpoint NO_(x) ^(sp) is calculated as a function of the reduction factor ξ.
 28. A method as claimed in claim 20, wherein: a NOx emission reduction factor ξ is set and the NOx emission setpoint NO_(x) ^(sp) is calculated as a function of the reduction factor ξ.
 29. A method as claimed in claim 13, wherein: determination of the dynamic torque setpoint of the diesel engine is obtained by inverting the NOx emission model.
 30. A method as claimed in claim 14, wherein: determination of the dynamic torque setpoint of the diesel engine is obtained by inverting the NOx emission model.
 31. A method as claimed in claim 15, wherein: determination of the dynamic torque setpoint of the diesel engine is obtained by inverting the NOx emission model.
 32. A method as claimed in claim 17, wherein: determination of the dynamic torque setpoint of the diesel engine is obtained by inverting the NOx emission model.
 33. A method as claimed in claim 20, wherein: determination of the dynamic torque setpoint of the diesel engine is obtained by inverting the NOx emission model.
 34. A method as claimed in claim 24, wherein: determination of the dynamic torque setpoint of the diesel engine is obtained by inverting the NOx emission model.
 35. A method as claimed in claim 13, wherein: the static torque setpoint T_(eng,s) ^(sp) of the diesel engine is obtained from a drive system torque setpoint T_(pwt) ^(sp) by an energy supervisor.
 36. A computer program product downloadable from a communication network and/or recorded on a computer readable medium, comprising program code instructions for implementing a method of controlling a hybrid drive system including at least one electric machine and at least one diesel engine wherein a static torque setpoint T_(pwt) ^(sp) of the diesel engine is acquired when the program is executed on a controller comprising: constructing a nitrogen oxide (NOx) emission model allowing estimation of the emissions of the diesel engine as a function of a torque T_(eng) of the diesel engine; calculating a dynamic NOx emission setpoint NO_(x) ^(sp) from the NOx emission model and from the static torque setpoint T_(eng,s) ^(sp) of the diesel engine; determining a dynamic torque setpoint of the diesel engine T_(eng) ^(sp) from NOx emission setpoint NO_(x) ^(sp); determining a dynamic torque setpoint of the electric machine T_(mot) ^(sp) as a function of the dynamic torque setpoint of the diesel engine T_(eng) ^(sp); and applying the dynamic torque setpoints to the drive system.
 37. A hybrid drive system for a vehicle comprising: at least one electric machine, at least one diesel engine and a controller, the controller executing a program performing a method comprising: constructing a nitrogen oxides (NOx) emission model allowing estimation of the emissions of the diesel engine as a function of a torque T_(eng) of the diesel engine; calculating a dynamic NOx emission setpoint NO_(x) ^(sp) from the NOx emission model and from the static torque setpoint T_(eng,s) ^(sp) of the diesel engine; determining a dynamic torque setpoint of the diesel engine T_(eng) ^(sp) from NOx emission setpoint NO_(x) ^(sp); determining a dynamic torque setpoint of the electric machine T_(mot) ^(sp) as a function of the dynamic torque setpoint of the diesel engine T_(eng) ^(sp); and applying the dynamic torque setpoints to the drive system. 